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TECHNICAL MEMORANDUM 


A TREETOPS SIMULATION OF THE STABLE MICROGRAVITY 
VIBRATION ISOLATION SYSTEM 

I. INTRODUCTION 


As a research facility for microgravity (pg) science, the International Space Station (ISS) will be 
used for numerous investigations such as protein crystal growth, combustion, and fluid mechanics 
experiments which require a quiescent acceleration environment across a broad spectrum of frequencies. 
Examples of the acceleration requirements for these processes are shown in figure 1 . Shown are the 
maximum magnitudes of acceptable accelerations, measured relative to nominal acceleration of gravity 
at the Earth’s surface, versus frequency, and the expected acceleration environment in the experiment. 
Note that these experiments are most sensitive to low frequency accelerations and can tolerate much 
higher accelerations at higher frequency. However, the anticipated acceleration environment on the ISS 
significantly exceeds the requirements shown in figure 1 . The ubiquity of vibratory acceleration distur- 
bance sources and the difficulty in characterization of these sources precludes source isolation, requiring 
vibration isolation to attenuate the anticipated disturbances to an acceptable level in the experiment. 



Estimate of Typical 
Space Station 
Environment 


Space Station 
Design Requirement 


Figure I . Microgravity acceleration requirements. 


The primary sources of vibration on the ISS can be cttegorized into three characteristic fre- 
quency ranges. At low frequencies. ~10~ 3 Hz, the dominant ; ccelerations are caused by gravity gradients 
and atmospheric drag. At high frequencies, above ~1 Hz, the vibrations are caused by sinusoidal steady- 
state sources such as pumps, compressors, electric motors, and fans, as well as transient sources such as 
impacts, astronaut motion, and high frequency components of thruster firings. The third characteristic 
frequency range of vibrations is the intermediate range of -10-3 to 1 Hz. The sources of accelerations in 
this range are mostly transient in nature, such as the motion t f astronauts and payloads around the ISS, 
as well as motion of the ISS caused by thrusters. Because of (heir transient nature, the effect of these 
vibrations on many experiments is difficult to analyze. The calculation of the resultant accelerations of 
the ISS are also complicated by the interaction of these vibration sources with the structural modes of 
the ISS, at least at the upper end of this frequency range. In the high frequency range, passive isolation 
techniques are often adequate to provide sufficient attenuation of vibration disturbances. However, 
isolation of low and intermediate frequency vibrations is not possible with passive isolation and requires 
active isolation. Hence, the development of active isolation systems is necessary to provide a quiescent 
acceleration environment as required by many microgravity science investigations. 

The basic objective of a vibration isolation system is to constrain the acceleration environment in 
the experiment by attenuating accelerations transmitted from the base to the isolated platform via um- 
bilicals or other direct disturbances. Base motion may be due to several sources such as crew motion, 
vehicle attitude control, or mechanical systems. In addition, disturbance forces which are transmitted 
directly to the platform (independent of the umbilicals) may result from crew contact or payload- 
generated sources such as pumps, fans, motors, and structural vibration of the isolated experiment. The 
required attenuation can be derived from the anticipated disturbance environment and required accelera- 
tion levels as shown in figure l. To provide the desired environment requires that the isolation system 
pass through the quasi-steady accelerations while providing attenuation above 0.01 Hz. At frequencies 
above 10 Hz, the required attenuation level is -60 dB, or 3 orders of magnitude of attenuation. To 
accomplish this, isolation in the presence of realistic umbilicals while rejecting direct disturbances 
requires an active isolation system. By sensing relative position and absolute acceleration of the plat- 
form, the control system causes the platform to follow the very low frequency motion of the base while 
attenuating the base motion above 0.01 Hz. High bandwidth acceleration feedback effectively increases 
the payload mass for disturbance rejection. Demonstration of :his level of performance in 6 degrees of 
freedom (DOF) using hardware cannot be accomplished on the ground due to gravitational coupling, but 
requires testing in a microgravity environment. Long periods >f experimentation are necessary to char- 
acterize the low frequency behavior, which is the most critical frequency range for active vibration 
isolation. 

For several years, an extensive effort has been underw ly in the aerospace community to develop 
active isolation systems for microgravity science experiments A joint effort by the Marshall Space 
Flight Center (MSFC) and McDonnell Douglas Aerospace (MDA) culminated in the design and flight 
test of the Suppression of Transient Accelerations By Levitation (STABLE) system. STABLE was the 
first successful flight test of an active isolation device for miciogravity science payloads and was flown 



on STS-73/USML-2 in October 1995. 1 Aside from being the first active microgravity vibration isola- 
tion system to be successfully implemented in orbit, STABLE was also significant in its development 
process. In less than 5 months, STABLE was designed, fabricated, tested, qualified, and delivered to the 
Kennedy Space Center. Nor is it trivial that the process was jointly conducted by NASA/MSFC in 
Huntsville, AL, and MDA in Huntington Beach, CA. The accelerated schedule and distributed location 
of the design team made imperative the use of a high-fidelity simulation for analysis and design. 

In order to verify the control system design and estimate on-orbit performance of STABLE, a 
high-fidelity, nonlinear, multibody simulation was developed using TREETOPS. 2 TREETOPS is a time 
history simulation developed for analysis of the dynamics and control-related issues of multibody 
structural systems. The name ‘TREETOPS” is indicative of the tree topology of linked multiple bodies, 
each of which may be rigid or flexible, with translations and large angle rotations between each body. 
Kane’s method is employed in the derivation of the equations of motion, which are numerically inte- 
grated to generate the time history response of the system. Extensive control system modeling capabili- 
ties are incorporated in TREETOPS, including a host of sensors and actuators along with controller 
models in the form of block diagram (transfer function), state space (multi-input/multi-output), and user- 
defined continuous or discrete controllers. An interactive setup program allows a convenient, easy-to-use 
interface with the simulation for model definition, input data editing, and error checking. For more 
detailed information on the analytical formulation and modeling aspects of TREETOPS, the reader is 
referred to the user’s guide. 2 

This report is written to document the development of the TREETOPS simulation of STABLE 
used in support of control system design. Because of the aggressive design, development, and delivery 
schedule for STABLE, sufficient time was not available for development of an accurate control design 
model. Instead, each of the six independent control channels were designed, based on a 1-DOF approxi- 
mation, and then analyzed and refined using the high-fidelity STABLE TREETOPS model. Hence, the 
STABLE TREETOPS model was essential for the successful design and flight test of STABLE. Section 
II describes the STABLE flight hardware and control system and section III is a parallel description of 
the simulation implementation of STABLE. Simulation results and results of the performance analysis of 
the STABLE system using TREETOPS are presented in section IV with conclusions in section V. 
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II. STABLE SYSTEM DESCRIPTION 


A. STABLE Hardware 

STABLE is comprised of a middeck locker, an isolated platform on which a thermal fluid con- 
vection experiment is mounted, three actuator assemblies, nine acceleration sensors, three position 
sensors, and associated electronics and control boards. The major elements of STABLE are shown in 
figure 2 (taken from ref. 1 ). The platform is suspended from the base of the locker box by three electro- 
magnetic actuators developed by MDA. The only physical ccnnections between the isolated platform 
and the base are flexible umbilicals that provide power and d ita to and from the platform. 



Figure 2. Configuration of STABLE system. 

Each actuator assembly provides two axes of force with a gap that allows ±1 cm of travel in each 
axis. This “sway space" is required to provide isolation from the very low frequency motion of the 
vehicle. The actuator is comprised of a paddle-shaped armatuie mounted to the base of the locker and a 
horseshoe magnet composed of an iron horseshoe and two Samarium-Cobalt pole pieces that produce a 
concentrated magnetic field in the gap between the poles. The armature paddle is inserted between the 
two permanent magnets of the horseshoe bracket as illustrated in figure 2. The coil assembly contains 
two independent coils, oriented such that when a current is passed through them, two independent 
orthogonal forces are produced in the plane, normal to the fiel i lines of the permanent magnet. By 
varying the current in the armature coils, a force is induced on the platform that is used to attenuate 
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disturbances and provide a quiescent acceleration environment on the platform. A high frequency accel- 
eration feedback control loop and a low frequency position feedback control loop are implemented to 
produce the required control force. The control algorithms are described in the next section. 

Six accelerometers and three relative position sensors are used to sense the isolated platform 
motion. The accelerometers are mounted in pairs on each of the three actuator mounting brackets and 
oriented to measure acceleration along the actuator force directions. Three additional accelerometers are 
mounted on the base and provide a measure of the base environment. The accelerometers used on both 
the platform and the base are model QA-2000, made by Sunstrand. Each of the three position sensors 
measures relative position of the platform with respect to the base at its particular location in two or- 
thogonal axes. It is composed of a laser illuminator, mounted on the platform, that projects a laser image 
across the sway space onto a photoresistive detector fixed to the base. The position of the platform with 
respect to the base is determined by processing the signals from the detector which indicate the location 
of the laser image on the face of the detector. 

B. STABLE Control Algorithms 

The key to the robust performance of STABLE is its six independent position and acceleration 
loops which provide high bandwidth acceleration feedback along with a positioning system that is 
insensitive to drift. A block diagram of this system is shown in figure 3. 



Figure 3. STABLE general block diagram. 
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The low bandwidth digital position controller uses measurements from the position sensors to 
compute the 6-DOF displacement of the floating platform and keep the floating platform centered in the 
sway space over a period of minutes. The acceleration loop is an analog controller with an ~50-Hz 
bandwidth to null the sensed acceleration of the platform. 

1. Position Controller 

A vibration isolator must attenuate “high-frequency” vibrations and be able to move with respect 
to the support structure and, thus, maintain an inertial position (or velocity) while the surrounding 
structure is in motion. To accomplish this, space must be provided around the isolated structure for it to 
“sway back and forth. The geometry of the sway space determines the lower frequency limit for attenu- 
ation of base motion. Below this low frequency limit, quasi-s-eady forces must be transmitted to the 
platform so that the platform will follow the low frequency motion of the support structure. The position 
controller serves this purpose. 

The position loop is a digital proportional-integral-derivative (PID) controller which uses sensor 
measurements to compute the position of the actuator with respect to a nominal centered position. 
Integral control is used to compensate for the error in calibration of accelerometer bias, unknown um- 
bilical bias force, and accelerometer bias drift due to temperature variation. Since the position sensors 
are not coincident with the actuators, they do not provide a di ect measure of displacement in the input 
axis of the actuator, and additional computations are performed in the digital processor to determine the 
position errors at the actuator gap. The algorithm for computing relative position at the gap from remote 
position sensor measurements is described in appendix A. Once computed, they are passed to the posi- 
tion controller, also in the digital processor, which calculates ; cceleration commands to each actuator. 
These acceleration commands are summed with the accelerometer signals and form the error signal 
input to the acceleration loop control law. 

The position control law operates in either of two modes, depending on the calculated actuator 
gap. While the gap position error is small, low PID gains are used to provide a very small restoring force 
which will not violate the acceleration requirements. If the gap error becomes large (meaning the plat- 
form is nearing the boundaries of its sway space), a set of high gains effectively increases the spring 
constant of the control law to prevent the platform from makii g mechanical contact. A block diagram of 
the position control loop is given in figure 4. 
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Coarse Bias 

Figure 4. Position loop controller block diagram. 


Figure 5 illustrates the logic employed for mode selection. Given the set of position error 
measurements, the high or low gain mode is determined by checking the magnitude of each error 
Sj (/' = l,---,6) for all six channels, as follows: 

• If during the previous sample period the high gain mode was in effect and the magnitude 
of any <5, is less than radius-low (RL), switch to low gain mode. 

• If during the previous sample period the low gain mode was in effect and the magnitude 
of any Sj is greater than radius-high (RH), switch to high gain mode. 



Edge of Gap 

Figure 5. Mode switch regions for position controller. 
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The gap error signal is input to the position controller which is passed through proportional, 
derivative, and integral channels. Two digital-to-analog (D/A) converters with different resolutions are 
used for the output. Fine resolution output signals are sent through a 12-bit D/A converter with a resolu- 
tion of 0. 1 |ig (200 jig full-scale) and coarse output signals use a larger range 1 2-bit D/A converter with 
a 10 jig resolution (20,000 (ig full-scale). The derivative charnel is low pass filtered to reduce high- 
frequency noise generated by differentiation. Logic is emplojed in the integral channel to reduce the 
excess that may be present due to buildup of large errors. While in low gain mode, the integral channel 
signal may become larger than can be adequately dealt with by the fine resolution output D/A converter, 
and thus the channel may be incrementally unloaded through the coarse channel D/A converter. The 
integral channel signal is compared with a maximum allowable value and if the signal is too large, an 
integer multiple of the coarse D/A converter least significant bit ( 1 0 jug) is subtracted from the fine 
channel and added to the coarse channel output. The integral channel for both fine and coarse outputs 
are checked for saturation prior to summing with other signals. 

2. Acceleration Controller 


The analog acceleration controller attempts to mitigate acceleration disturbances sensed by the 
accelerometers. The analog accelerometer measurement forms the acceleration error signal input to the 
analog acceleration controller. A rate feedback-type acceleration controller provides stability robustness 
since the actuator and sensor are spatially collocated. A low pass filter provides rolloff at a nominal 
bandwidth of 50 Hz. Performance of the acceleration loop is limited by controller bandwidth, acceler- 
ometer noise, resolution- and temperature -dependent bias variations, and disturbances transmitted 
through the umbilical connections. The output from the acceleration loop is the current command to the 
actuator coils. The acceleration control loop is implemented in TREETOPS with the transfer function 

. . . _ , K a (s + I0*2n) 

Acceleration Control = ~ j 0s( " + 2 ^) ~ < N / m / s2 ) * 

where K 4 = 5.5 x I0 4 . 

The next section describes the TREETOPS implementation of the STABLE dynamics and control 
algorithms. 
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III. DESCRIPTION OF STABLE TREETOPS MODEL 


In order to verify the control system and estimate the performance of STABLE in the orbital 
environment, a detailed structural and control model of STABLE was developed for TREETOPS. Since 
the platform of STABLE is floating freely within the STABLE locker box that is rigidly attached to the 
space shuttle, the STABLE system was modeled as two rigid bodies (first body for the space shuttle and 
rigidly attached locker box and second body for the isolated platform). The only physical connections 
between the platform and base are the umbilical cables (that transfer data and power between the plat- 
form and base) which are modeled as a 6-DOF hinge connection with spring stiffness specified to match 
the dynamic properties of the umbilical cables. 

The TREETOPS model of the STABLE system must accurately model all of the actual flight 
hardware components to precisely represent the flight system. Note that TREETOPS terminology is 
indicated by spelling with capital letters. The STABLE actuator that yields two orthogonal forces can be 
modeled using two TREETOPS JET actuators and each STABLE accelerometer sensor can be modeled 
using the TREETOPS ACCELEROMETER sensor. Actuator and sensor models are idealized and do not 
include high-frequency dynamics. Since there is no built-in TREETOPS position sensor model that 
exactly corresponds to the STABLE position sensor, a mathematical position sensor model was devel- 
oped using the built-in TREETOPS POSITION VECTOR sensor and STAR sensor. This position sensor 
model was implemented in a user-defined discrete controller subroutine and incorporated with the main 
TREETOPS dynamics simulation. The convention of the local coordinates and locations of the actuators 
and sensors are shown in figure 6. In this figure, Fjj, A,y, and e,y denote /th force component of /th actua- 
tor, /th acceleration component of /th accelerometer, and /th position error component of /th position 
sensor, respectively ( 1 st, 2nd, and 3rd components stand for .v, v, and r axis components, respectively). 



Figure 6. STABLE coordinate systems. 
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With supplied mass properties (mass and moments of inertia) of the two bodies and the locations 
of center of mass (CM), actuators, accelerometers, position srnsors, and umbilical connections, 
TREETOPS determines the kinematics and dynamics of the STABLE system. The fast, continuous- 
acceleration control loop was implemented in the STABLE TREETOPS model using the built-in TREE- 
TOPS continuous block diagram controller and the slow digi al position control loop was implemented 
in the STABLE TREETOPS model using a user-defined disc ete controller subroutine. 

A. STABLE Structural Model 

The STABLE system, including the space shuttle, was modeled as two rigid bodies connected 
by linear and nonlinear quadratic springs. Since the objective of the STABLE controller is to achieve 
attenuation of an acceleration disturbance from the locker bo.; to the STABLE platform, the space 
shuttle and attached locker box was modeled as one arbitrary rigid body that gives a disturbance to the 
STABLE platform through umbilical cables. 

According to the TREETOPS tree topology, the space shuttle and attached locker box is defined 
by body #1 and linked by hinge #1 with 6 DOF (3 rotational and 3 translational) with respect to the 
origin of the inertial coordinate system. The platform floating inside the locker box is defined by body 
#2 and connected to body #1 through hinge #2 with 6 DOF. The umbilical connection between the 
platform and locker box is modeled as the combination of six linear spring devices and two quadratic 
spring devices with a 10-m undeformed length. This undefomed length was used to produce the desired 
linear force without introducing artificial rotational coupling. 

For body #1,12 nodal points are chosen to represent the CM, origin of local coordinate system 
of body #1, three position sensor detector attach points, the C vl of body #2, and six umbilical cord 
attach points. For body #2, 16 nodal points are chosen to represent the CM, origin of local coordinate 
system of body #2, three actuator attach points, six accelerom iter attach points, three position sensor 
laser attach points, and two umbilical cord attach points. Table 1 summarizes the nodes of body #1 (for 
example. B I N2 denotes node #2 of body # 1 ) and the nodes of body #2 are summarized in table 2 (for 
example, B2N2 denotes node #2 of body #2). The definitions 3 f all hinges of STABLE TREETOPS 
model are summarized in table 3. 
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Table I. Node definitions of TREETOPS STABLE body #1. 


Node 

Description 

Location in Body Coordinates (meter) 

B1N1 

CM of body #1 

(0,0,0.02) 

B1N2 

Origin of body #1 coordinate 

(0,0,0) 

B1N3 

Position sensor laser #1 

(0.1537,0,0) 

B1N4 

Position sensor laser #2 

(-0.1537,0,0) 

B1N5 

Position sensor laser #3 

(0 -0.1537,0) 

B1N6 

CM of body #2 

(0.004,-0.02,0.067) 

B1N7 

X umbilical #1 

(10.0686,-0.0787,-0.0205) 

B1N8 

/umbilical #1 

(0.0686,9.9213,-0.0205) 

B1N9 

Z umbilical #1 

(0.0686,-0.0787,9.9795) 

B1N10 

X umbilical #2 

(9.9314,-0.0787,-0.0205) 

B1N1 1 

/umbilical #2 

(-0.0686,9.9213,-0.0205 

B1N12 

Z umbilical #2 

(-0.0686,-0.0787,9.9795) 


Table 2. Node definitions of TREETOPS STABLE body #2. 


Node 

Description 

Location in Body Coordinates (meter) 

B2N1 

CM of body #2 

(0.004,-0.020,0.067) 

B2N2 

Origin of body #2 coordinate 

(0,0,0) 

B2N3 

Position detector #1 

(0.139,0,0) 

B2N4 

Accelerometer #1 

(0.0411,0.0942,0.0754) 

B2N5 

Actuator #1 

(0,0.1433,0.0423) 

B2N6 

Position detector #2 

(-0.139,0,0) 

B2N7 

Actuator #2 

(-0.1218,-0.1735,0.0423) 

B2N8 

Accelerometer #3 

(-0.0582,-0.1696,0.0754) 

B2N9 

Position detector #3 

(0,-0.139,0) 

B2N10 

Accelerometer #5 

(0.0582,-0.1696,0.0754) 

B2N11 

Actuator #3 

(0.1218,-0.1735,0.0423) 

B2N12 

Umbilical #1 

(0.0686,-0.0787-0.0205) 

B2N13 

Umbilical #2 

(-0.0686,-0.0787-0.0205) 

B2N14 

Accelerometer #2 

(0,0.0942,-0.0206) 

B2N15 

Accelerometer #4 

(-0.0918,-0.1388,-0.0206) 

B2N16 

Accelerometer #6 

(0.0918,-0.1388,-0.0206) 


Table 3. Hinge definitions of STABLE TREETOPS model. 


Hinge 

Connecting Nodes 

No. of DOF 

L1_in-L1_out 

L3_in-L3_out 

1 

B0N0-B1N2 

3RD0F, 3TD0F 

(1,0,0)-(1,0,0) 

(0,0,1)-(0,0,1) 

2 

B1N6-B2N1 

3RD0F, 3TD0F 

(1,0,0)-(1.0,0) 

(0,0,1)-(0,0,1) 






B. STABLE Sensors Model 


STABLE has six QA-2000 accelerometers on the Hotting platform to measure acceleration at 
the attached nodes. This accelerometer was modeled as TREETOPS ACCELEROMETER (AC) sensor 
that measures the acceleration of body at the specified node in the specified direction. Therefore, six AC 
sensors are attached on nodes #4, #8, #10, #14, #15, and #16 of body #2. STABLE has three position 
sensor assemblies to measure the relative position errors between the floating platform and locker base 
at three position sensor locations. The measured position errors are sent to the digital computer to com- 
pute the position error at the actuator gap using the derivation in appendix A. The relative displacement 
at the gap is the input to the position control law. Each STABLE position sensor assembly consists of a 
laser attached on the platform and an image detector attached on the locker base. The image detector 
measures two orthogonal displacements projected on the detector by the laser. 

Since there is no built-in TREETOPS sensor corresponding to the STABLE position sensor 
assembly, a mathematical model of a corresponding position ; ensor assembly was developed using the 
TREETOPS POSITION VECTOR (P3) sensor and STAR (ST) sensor. For example, to model the 
STABLE position sensor assembly at #1 position, a P3 sensor is added on the laser location of the 
platform to measure three displacement components between the laser attached on the platform and the 
image detector attached on the locker base, i.e., <y p\ p \ ). Two ST sensors are also added at the 
laser location of the platform and the detector location of the locker base to measure two relative rota- 
tional angles of the laser with respect to the detector about r aid y axis, i.e., (0\ - ,Q \ V ). Then, the output 
of the TREETOPS position sensor #1 corresponding to the STABLE position sensor assembly #1 can be 
given by 


• v i - -y P i + x P \ 


*ta/U) 


1 _- 


■p 1 


l P 1 


*tan6 „ . 


Similarly, the outputs of the TREETOPS STABLE position seisor #2 and #3 can be given by 

V ’2 = y pi - Xp 2 1 * tan 6 2-_ - 


-2 - pi 


X P 2 


*tandi Y , 


3 ~ ~ x p3 + p3 * tan Oj 


-3 - - p 3 


|>>3 


* ran 6 


t.Y 


0) 

( 2 ) 


(3) 

(4) 

( 5 ) 

( 6 ) 
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Although not part of the STABLE flight hardware, for the purpose of performance evaluation, 
one P3 sensor is attached on the CM of body #1 (B 1 N1 ) to measure the relative movement of the CM of 
the platform with respect to the locker base. Also, two IMU SENSOR (IM) sensors are used to measure 
the Euler angles of local frames of bodies #1 and #2. These sensors are summarized in table 4. 


Table 4. Definition of TREETOPS STABLE sensors model. 


Sensor ID No. (Type) 

Measurement Quantity 

Location (Direction) 

SE 101 (AC) 

Acceleration #1 

B2N4 (1,0,0) 

SE 102 (AC) 

Acceleration #2 

B2N14 (0,0,1) 

SE 201 (AC) 

Acceleration #3 

B2N8 (-0707,0.707,0) 

SE 202 (AC) 

Acceleration #4 

B2N15 (0,0,1) 

SE 301 (AC) 

Acceleration #5 

B2N10 (-0.707,-0.707,0) 

SE 302 (AC) 

Acceleration #6 

B2N16 (0,0,1) 

SE 401 (ST) 

e 2 &0, of laser #1 

B2N3 (0,1,0) & (0,0,1) 

SE 402 (ST) 

0 2 &e 3 of detector #1 

B1N3 (0,1,0) & (0,0,1) 

SE 501 (ST) 

0 2 & 0} of laser #2 

B2N6 (0,1,0) & (0,0,1) 

SE 502 (ST) 

0 2 &0 3 of detector #2 

B1N4 (0,1,0) & (0,0,1) 

SE 601 (ST) 

0 ! & 0 3 of laser #3 

B2N6 (1,0,0) & (0,0,1) 

SE 602 (ST) 

0i &0 3 of detector #3 

B1N4 (1,0,0) & (0,0,1) 

SE 701 (P3) 

Xf)i-Yp\ -Zp\ 

(B1N3)-(B2N3) 

SE 801 (P3) 

Xpl'Ypl'Zpi 

(B1N4)-(B2N6) 

SE 901 (P3) 

XpyYpyZpi 

(B1N5)-(B2N9) 

SE 91 1 (P3) 

Relative movement of body #2 
CM w.r.t body #1 

(B1N6)-(B2N1) 

SE 921 (P3) 

Movement of body #1 CM 

(B1N6)-(B2N1) 

SE 912 (IM) 

Euler angles of body #2 frame 

B2N1 

SE 922 (IM) 

Euler angles of body #1 frame 

B1N2 


C. STABLE Actuator Model 

The STABLE has three MDA dual-axis actuators which generate two orthogonal forces and are 
modeled as two TREETOPS JET (J) actuators. The actuator applies a force at the node of the platform 
where the actuator is attached. Each JET actuator gives constant force that is defined using a FUNC- 
TION GENERATOR (FU) at the specified node to the specified direction. Therefore, six JET actuators 
are attached on nodes #5, #7, and #1 1 of body #2. For the purpose of performance analysis, a distur- 
bance can be given by applying a force on the CM of body #1 using a separate JET actuator. These 
actuators are summarized in table 5. 
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Table 5. Definition of TREETOPS STABLE actuators model. 


Actuator ID No. (Type) 

Measurement Quantity 

Location (Direction) 

AC 101 (J) 

Actuator force #1 

B2N5 (1,0,0) 

AC 102 (J) 

Actuator force #2 

B2N5 (0,0,1) 

AC 201 (J) 

Actuator force #3 

B2N7 (-0.707,0,707,0) 

AC 202 (J) 

Actuator force #4 

B2N7 (0,0,1) 

AC 301 (J) 

Actuator force #5 

B2N11 (-0.707,-0.707,0) 

AC 302 (J) 

Actuator force #6 

B2N11 (0,0,1) 


D. STABLE Controller Implementation 

The high performance characteristics of the isolator are the result of active feedback loops 
involving actuators, sensors, and electronics. The isolated platform is controlled by means of six 
independent control channels, one for each actuator force direction. Each STABLE control channel 
consists of one fast inner acceleration control loop and one slow outer position control loop, as described 
in section II. 

In order to implement this STABLE controller in the TREETOPS dynamics model, three 
TREETOPS control modules (CBDC, USDC, and USCC) were used. The inner acceleration controller 
loop was modeled using continuous block diagram controller (CBDC) and the outer position controller 
loop was modeled using user-defined discrete controller (USDC) subroutine. The user-defined continu- 
ous controller (USCC) subroutine was used to connect the inner acceleration controller loop and the 
outer position controller loop. The USCC and USDC subrouti les are documented in reference 3. 

1. Position Control Loop 

The STABLE position PID controller was implemented in the TREETOPS STABLE model 
using a USDC subroutine. ^ This USDC routine generates the ; cceleration command, based on position 
error, and prints out the intermediate results for the purpose of performance analysis. This routine also 
provides a procedure that determines the output of the TREETOPS position sensor model corresponding 
to the STABLE position sensor assembly using equations ( 1 ) — 6). The relative displacements at the 
actuator gap are computed from the position sensor measurements in a two-step process. First the dis- 
placements of the platform CM —eg '^x <x •@y ct , ) are determined and then the displace- 

ments at the actuator gaps are determined, given the location of the CM, the sensors, and the actuators. 
The relationships between the outputs of STABLE position sensors, the relative displacements of the 
platform CM, and the relative displacements at the actuator gaes are derived in appendix A. The discrete 
position control loop computes an acceleration command which will force the position error at the 
actuator gap to tend to be zero. 
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2. Acceleration Control Loop 


Prior to implementing the acceleration control law in CBDC, the USCC sets the control gain and 
performs part of the acceleration control law computation. Inputs to the USCC are the high or low gain 
flag (determined by the mode switching logic in the position control law), acceleration command from 
the position loop ( V f ), and the measured acceleration (.¥)• USCC is used to select and output the proper 
(high or low) acceleration control gain ( K A ) and products ( K A * x) and ( K A *x c ). These quantities are 
output for use by CBDC. (In the flight system, the acceleration command and measurement are summed 
to form the error signal which is multiplied by the acceleration gain and processed by the analog accel- 
eration control law. To realize this operation in TREETOPS controller block structure, the gain is multi- 
plied by acceleration command and measurement prior to summation.) 

The STABLE inner acceleration controller has six independent channels for six actuator forces 
and is modeled using TREETOPS CBDC. STABLE employs an analog single-input/single-output 
acceleration control architecture. Based on the collocated acceleration measurement, a control force is 
computed to reduce the inertial motion of the platform. Each channel has the same control law and 
continuously determines control force to compensate for an acceleration error at the actuator attach 
point. CBDC generates a control force for each actuator force direction by summing the outputs of 
USCC (K a *x j and (K A *.v ( .)and processing the error signal by the transfer function representing the 
analog control circuitry. 

For the TREETOPS model, a total of thirty inputs, five for each channel, are read into CBDC. 
The inputs of each channel are two FUNCTION GENERATORS (FU) that provide a dc bias and white 
noise acceleration component, an acceleration measured by the AC sensor, and two outputs of USCC, 

K a *.v and K A *x c . CBDC generates a control force for each actuator output direction by computing 
(-K A *x-K a *x c ) and processing the error signal by the transfer function representing the analog 
control circuitry. This CBDC generates twelve outputs, two for each channel — the actuator control force 
and the total acceleration from acceleration control loop. 
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IV. SIMULATION RESULTS 


This subsection describes numerical results of STABLE stability and attenuation performance 
analysis obtained from TREETOPS simulation. The input file of the STABLE TREETOPS simulation is 
documented in reference 3. For the stability analysis, two cast's are considered: One is to check stability 
of the STABLE system by measuring the transient response of the platform with given initial movement 
of platform and the other is to check stability of the STABLE system by measuring the transient re- 
sponse of the platform with a given pulse-type disturbance on base. The attenuation performance analy- 
sis was done to determine the acceleration attenuation curve by measuring the acceleration at the plat- 
form CM due to the various frequencies of sinusoidal disturbances given to the base. Mass properties of 
the STABLE platform are described in table 6. 


Table 6. Mass properties of the STABLE platform. 


Mass (Kg) 

12.646 

A.v A-v At- At. A.< (Kg-m-i 

0.167,0 108, 0.179, 0, 0,0 


As explained in section III, the umbilicals are modeled as combinations of linear and quadratic 
spring devices. Based on the measured spring stiffness, the linear spring stiffness coefficients used for 
TREETOPS STABLE model are 18 N/m in the positive X-axis direction, 13.5 N/m in the positive F-axis 
direction, and 20 N/m in the positive Z-axis direction. The quadratic spring stiffness coefficient used for 
TREETOPS STABLE model is -672 N/m 2 to the positive Z-a; is direction. 


Final control parameters used for the acceleration conti oiler and PID position controller of the 
flight STABLE system were determined through iterative design and performance simulation and are 
summarized in table 7. In figure 4, K tL and K m correspond to he low and high gain mode values for X/, 
respectively. Similarly, K and correspond to the low and high gain mode values for K respec- 
tively. RH and RL are parameters which specify the mode switch regions in figure 5. 
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Table 7. STABLE control parameters. 


Control Parameters 

Low Gain Mode 

High Gain Mode 

K A (N*se& tn) 


5.5e+4 

Lead-lag filter 

1 

1 f s + 2k *10 ^ 

Transfer function 

lol .v+27t ) 

lol ,v + 2k J 

Kd 1/sec) 

8.8e-2 

2 

Kp( 1/sec 2 ) 

4e-3 

4e— 1 

/(/(I /sec 3 ) 

4e-5 

4e-2 

Max. Cuing) 

20 

- 

RL (mm) 

2 

- 

RH(mm) 

- 

8 

LSB (|ig) 

0.2 

10 

AC,/(MQ) 

10 

- 

«l 

2 

1 

Ts 

0.1 

0.1 

N S 

- 

100 


A. Case 1: Transient Response Analysis With Initial Displacement 

In order to investigate the stability of the STABLE control system with acceleration and position 
control, the STABLE platform was initially displaced from the nominal resting position by 10 mm in 
each axis direction. Accelerometer bias and white measurement noise are added to the accelerometer 
outputs to represent the hardware characteristics of the STABLE accelerometers. High gain mode was 
set initially and the transient response was then calculated from the TREETOPS simulation. The accel- 
eration biases for the six accelerometers were chosen arbitrarily for this simulation, as shown in table 8. 


Table 8. Acceleration biases of six accelerometers. 


Acc. Bias 

All 

A13 

A21 

A23 

A31 

A33 


105 

-155 

85 

-125 

25 

115 


The white acceleration measurement noise was generated with random numbers multiplying the follow- 
ing transfer function implemented in CBDC: 

_ + a) 

M h a(s 2 + 2 %(Q n s + (Q n 2 ) ’ 

where K=2x\0~^ m/sec% a=2rci20) rad/sec, to„=2/r( 1 00) rad/sec. and £=0.85. 
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Numerical simulation results of the STABLE transient response analysis are shown in 
figures 7—14. The gains of the PID controller were selected based on the switching logic and are shown 
in table 7. 


Figure 7 shows the following time responses: The output of the six PID position controller’s 
channels; the accelerations A 1 1 , A2 1 , and A3 1 are parallel to the direction of accelerometers # 1 , #3, 
and #5, respectively; and the accelerations A 13, A23, and A33 are in the Z-axis direction parallel to the 
direction of accelerometers #2, #4, and #6, respectively. As shown in this figure, high gain mode was 
initially active for about 40 sec and when all six channels satisfied the low gain criteria for more than 
N s time samples simultaneously, the position controller changed into low gain mode to get precise 
control of position. It should be noted that all six channels show outputs with equal magnitude and 
opposite signs of the accelerometer biases given in table 8. 


Each channel of the position controller has two output signals — coarse acceleration (PC) quan- 
tized with 10 pig least significant bits (LSB) and fine accelerat on (PF) quantized with 0.2 pg LSB. 

In figure 8, the coarse and fine channels are designated PCij and PFij, respectively. Defining A\ and A-> 
be acceleration output without and with quantization, respectively, the equation of quantization may be 
given as follows: 



(\ „ 1 LSB^ 


int 

U, + 

USB 


V 1 2 J 



LSB * sign( A j ) 


( 8 ) 


The sum of PF and PC is the output of the PID position controller and is fed back to the acceleration 
controller implemented in CBDC at every sampling time. These results are shown in figures 8 and 9. 


Detailed outputs of the sixth channel of the position coitroller are presented in figure 10. In this 
figure, the first row of plots show position error, position error rate, and integrated position error in the 
Z-axis direction at the gap of actuator #3. Multiplying these error signals by the PID gains given in table 
7 yields cp, cr, cijow, and ci_high. The coarse acceleration output (PC) is the sum of cp, cr, and 
ci_high. The fine acceleration output (PF) is the sum of cp, cr, and cijow. The last row of plots shows 
indication of PID position high or low gain mode. 


Figure 1 1 shows the outputs of six accelerometers: A 1 1 , A 1 3, A2 1 , A23, A3 1 , and A33 are the 
measured accelerations of the accelerometers #1, #2, #3, #4, #3, and #6, respectively. It is noted that the 
acceleration controller immediately senses abrupt changes of several hundred microgravity accelerations 
and reduces the acceleration levels to the given output noise le/el of accelerometer. The acceleration 
command from the position controller results from the positior controller command to move the 
STABLE platform from the initially displaced position to the nominal rest position; therefore, the accel- 
eration controller is operating as designed. 


Figure 12 shows the position errors measured by the ST ABLE position sensors described in 
section II. In figure 12, LI and Zl, Y2 and Z2, and Y3 and Z3 aie Taxis and Z-axis directional position 
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errors of the STABLE position sensors #1 , #2, and #3, respectively. It is shown that the PID position 
controller brings down position error of the STABLE platform from initial 10 mm to zero mm in a 
timely manner. 

The controller position errors at the CM of the STABLE platform, shown in figure 13. are calcu- 
lated using equations ( 1 1 )-(2 1 ) in appendix A together with the position errors measured by position 
sensors. To check the validity of equations ( 1 1 )-(2 1 ), the position errors at the CM of the STABLE 
platform were obtained from TREETOPS' built-in position sensor at the CM of STABLE platform and 
hinge #l’s output, and plotted using the dotted line in figure 13. These results confirm the validity of 
equations (1 1H21) and the effectiveness of the STABLE PID position controller. 

Figure 14 shows the position errors at the actuators’ gap, determined using equation (28) in 
appendix A. In this figure, outer dotted lines denote RH (±8 mm) and inner dotted lines denote RL 
(±2 mm), and these results verify the switching logic for the high/low gains of the PID position 
controller explained in section II. 



Time (sec) 

Figure 7. PID position controller output with initial displacement. 
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Figure 9. Detailed PID position 
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Figure 10. Detailed PID position controller output of sixth channel with initial displacement. 
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B. Case 2: Transient Response Analysis With Initial Excitation 


In order to demonstrate how well the STABLE system can overcome sudden disturbance, a 
1,000-pg step disturbance input was given to the base of the Z-axis direction for 1 sec, and then a tran- 
sient response analysis was performed using input data tables 7 and 8 with initial low gain mode. The 
biased acceleration and white noise of the accelerometers used for this simulation are the same as those 
used in case 1 . 

The numerical results of this transient response analysis are shown in figures 15-22. Figures 
15-17 show the output of the six PID position controller’s channels with detailed output of sixth channel 
of position controller presented in figure 18. Figure 19 shows the measurements of the reference six 
accelerometers. As shown in these figures, the abrupt disturbance is sensed and the controller shifts from 
low gain mode into high gain mode immediately. The STABLE controller effectively responds to abrupt 
input accelerations and reduces the platform acceleration to the noise level of the accelerometers. It is 
also noted that the low gain mode may be recovered in about 40 sec. 

Figure 20 shows the position errors measured by the STABLE position sensors and the position 
errors at the CM of the STABLE platform are shown in figure 21. Figure 22 shows the position errors at 
the actuator gap. 


C. Case 3: Attenuation Performance Analysis Results 

The main objective of the STABLE system is to provide a low acceleration environment across a 
broad spectrum of frequencies using an active isolation controller. An attenuation performance analysis 
of the STABLE control system was performed and numerical results are presented in this section. The 
acceleration attenuation of STABLE system was estimated by measuring the accelerations at the plat- 
form CM for given sinusoidal disturbances of various magnitudes and frequencies. For this analysis, the 
accelerometer bias, accelerometer noise, and initial displacement are not included. The acceleration 
attenuation curve was determined by taking the ratio of acceleration magnitudes at the platform CM to 
the given disturbance acceleration magnitudes across the frequency range of 0.001 through 100 Hz. 
Figure 23 shows the calculated acceleration attenuation curve of the STABLE system. With a break 
frequency of approximately 0.01 Hz and rolloff of approximately 20 dB/decade, good attenuation is 
achieved. The slight amplification below 0.02 Hz is a consequence of low frequency acceleration com- 
mands produced by the position controller to keep the platform centered in the sway space. 
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Figure 16. Detailed PID position 
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Figure 21. Position error at the CM of STABLE platform with initial excitation. 
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V. CONCLUSIONS 


This report documents the details of the control/dynamics simulation used to design the control 
laws and define the expected performance of the STABLE isolation system. The multibody simulation 
tool, TREETOPS, was used to implement the 6 DOF dynamics of the floated platform as well as the 
dynamical models of all sensors, actuators, and the control algorithms in the software. Where possible, 
models were based on measured data and transfer functions; and, therefore, the simulation should be a 
good representation of the actual hardware performance. It should be noted that the performance predic- 
tion in terms of an attenuation transfer function does not meet the requirement that was established for 
STABLE. This was due to a situation in which, because of schedule pressure, it was necessary to use 
hardware that was designed for another program, hardware whose electronics characteristics forced a 
compromise in the design of the bandwidth of the position control loop. The results in this report will be 
compared to flight data and will be published in a separate report. 
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APPENDIX A — Derivation of Relative Position at the Actuator Gap 


Since the six control loops are decoupled and use measurements collocated with the actuators, 
the input to the position control law is the relative displacement of the isolated platform with respect to 
the base at the actuator gap and in the direction of actuator force axes. However, the position sensors are 
located remotely from the actuators in the STABLE configuration. A laser mounted on the isolated 
platform illuminates a detector fixed on the base to provide a measure of the relative motion. Each 
detector measures two axes of motion in the plane of the detector, normal to the incident laser beam. 
These six relative position measurements are used to compute the relative motion at the actuator gaps 
based on the locations of the actuators, lasers, and detectors through a two-step process. First, the 
6-DOF motion of the CM is computed from the sensor measurements and then the resulting displace- 
ment at the actuator gaps is determined. The equations defining these two steps are derived in the fol- 
lowing sections. 


A.l Center of Mass Relative Displacements 


Figure 24 indicates the vectors and coordinate frames used for the computation of the position 
error projected onto the detector plane of sensor #1 that results fron motion of the platform CM, 

( -'Vi,- ■ x 'ci> '-eg Ay, Ay, A \ The relative position, as measured by the detector at sensor #/ (/=1 ,2,3), 
is the projection of the error vector, e,, onto the detector plane. The error vector is given by 


£/ — ZI() + Can + P + £/ ~ (5 / + £ki ) » 


(1) 


where the T frame is fixed in the base such that the center of the detectors lie in the Y_\ fSplane; the Z 
frame is fixed in displaced platform such that the axes of the Z and Y frames are parallel and the origins 
coincident when the platform is in the null position; s' i locates the /th laser in the fixed Y frame; S ■ 
locates the /'th laser in the displaced Z frame; A, is the vector defin ng the location of the /th detector 
with respect to the /th laser in the null position; locates the CM n the Y frame; r cm is the displace- 
ment of the CM with respect to the Y frame; and p locates the origin of the displaced platform frame 
with respect to the CM in the Z frame. For the STABLE configurat on, position sensor locations are 
shown in figure 25, from which these vectors are defined as 



+ r 2^-2 + * 3—3 

.v Y i 

+ V Y i 4 - z 

Cg - 1 

'O'-- eg- 

~>]Z\ 

— /% Z 2 — / ^ 
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S | — 5 ] Z | ; S 2 = ~$2 — 1 ; 


S 3 =-S 3 Z 2 


S\=S l Y l ; S , 2=-S 2 Y l -, S'3=-S 3 Y 2 

A l =A l Y 1 ; A 2 =-A 2 Y\\ A 3 =-A 3 Y 2 . 



Displaced CM Location 


Figure 24. Definition of coordinate frames for relative position measurements. 
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Figure 25. Position sensor detector locations. 
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Assuming small angles, the transformation from the Z frame to the Y frame is 


T YZ ~ 


1 - 0 . 




e. 
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-e 


>v. 


x c* 


-e 


- V 0! 


e 




Realizing all vectors in the Y frame yields the error vectors 
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(5) 


Figure 26 shows the projection of onto sensor #1, for which the detector plane is the 
Yi K 3 plane. Note from figure 25 that the sensor frame is not aligned with the Y frame, but T\ j = —Y\ 
and fp = -Yi. The error vector originates at the center of the detector and terminates at the laser in the 
displaced platform frame. 
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This nonlinear equation can be solved for the platform CM motion, ( A cg , ycg’“c t (' , ^.v ( . ( ,,^v ) given 

the six sensor measurements using the following relationship: 


9. — — (£|2 + £?? ) / (2 * (s + <0) 


i'V 


(ID 


0 V — — (£] 3 — £ 33 ) / (2 * (.v + ■/)) 


( 12 ) 


r, = -£33 + (£ 13 +£ 23 )/ 2 + £32 * 9 V - r 3 * d Y 2 4 (53 + d + r 2 ) * 0 V * 6 . 


* 


* 


Of CV 


(13) 


C-) = (£-n — £1 7 ) / 2 + /■] * 9 , + 0. * (—£37 + / 3 * 6 — (s 3 + J + r 7 ) * ) 

*" 1 “ Tv Tv ‘ ■" - TV ' “ TV 


(14) 


A, =1+0. 
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(15) 


A~> — (.93 + d) * A] + £3 * 9 x * 0 . — C 7 


• <V ~i V 


(16) 
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a 3 = r 3 * r { - r 2 * (5 3 + d) 


(17) 


y cg = (~A 2 +i]A%- 4*A 1 *A 3 )/(2*A 1 ) 
x ce = ~ e M + r 3 * 6 x -{s 3 +d + r 2 )*0. -v *0- 

A - *r« “ct? “c» 

0 V =(r,+0 *0. *y )/{s 3 +d + y ) 

' ri* ■ Yt> Vtf tS ' A 

: ci? = < e l3 +e 23 )/2 + r ‘y *e -r,*0 -0 *0 

A ~ *- Ci> ■ (- ’ Cl? ■ Ci> ' 


where 5j - S 2 =s-, S 3 = sy, and A i - A 2 = A 3 =d . 


A.2 Relative Displacements at the Actuator Gap 


(18) 

(19) 

( 20 ) 
( 21 ) 


The displacements at the actuator gap may now be determined from the displacements of the 
platform CM computed using equations (1 1 )— (2 1 ). Figure 27 indicates the same coordinate frames and 
vectors locating the CM as figure 24 and also introduces the vector defining the location of the center of 
the gap associated with actuator#/ (/= 1,2.3), designated by R, in the Y frame and by R- in the Z frame. 
From figure 27, the relative displacement measured at the gap of actuator #z is 

&i = ~&i + Lq + Lem + P + /? / • (22) 
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Figure 28. Actuator coordinate frames. 
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Equation (22) may now be evaluated for the relative displacement at the gap of actuator #/ 
in the Z coordinate frame: 
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After simplifying. 
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( 27 ) 


The gap displacements must be expressed in the appropriate actuator frame using the transformation 
matrices defined in equation (26) and the immediately preceding text. In the actuator frames, 

-i' ~ T s . z &i - which results in 
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